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Abstract 

We discuss some generic properties of lepton number violating ) pro- 
cesses and their relation to different entries of the Majorana neutrino mass 
matrix. Present and near future experiments searching for these processes, 
except the neutrinoless double beta decay, are unable to probe light (eV 
mass region) and heavy(hundred GeV mass region) neutrinos. On the other 
hand due to the effect of a resonant enhancement, some of Ijj decays can 
be very sensitive to the intermediate mass neutrinos with typical masses 
in hundred MeV region. These neutrinos may appear as admixtures of 
the three active and an arbitrary number of sterile neutrino species. We 
analyze the experimental constraints on these massive neutrino states and 
discuss their possible cosmological and astrophysical implications. 
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1 Introduction 



Recent evidences for neutrino oscillations leave nearly no room for doubts that 
neutrinos are massive particles. After all this point of view is becoming con- 
ventional. Solar neutrino deficit, atmospheric neutrino anomaly and results of 
LSND neutrino oscillation experiment all can be explained in terms of neutrino 
oscillations implying non-zero neutrino masses and mixings [Q. On the theoreti- 
cal side there are also many indications in favour of non-zero neutrino masses 
following from almost all phenomenologically viable models of the physics beyond 
the standard model(SM). These models typically predict Majorana type neutrino 
masses suggesting that neutrinos are truly neutral particles. 

The neutrino oscillation searches fix both the neutrino mass square difference 
Smfj = mf—m'j and the neutrino mixing angles, leaving the overall mass scale and 
the CP-phases arbitrary. Since the latter has no effect on neutrino oscillations, 
the important question of whether neutrinos are Majorana or Dirac particles 
cannot be answered by these searches. 

Majorana masses violate total lepton number conservation by two units AL = 
2. Thus lepton number violating(^ ) processes represent a most appropriate 
tool to address the question of the Majorana nature of neutrinos. A celebrated 
example of ^ -process, most advanced experimentally and theoretically, is the 
neutrinoless nuclear double beta (Oi//3/3) decay (for a review see [@, H). The Oz//5/5- 
experiments achieved unprecedented sensitivity to the so called effective Majorana 
neutrino mass {m,y)ee 0, which in the presence of only light neutrinos coincides 
with the entry of the Majorana neutrino mass matrix {m,^)ee = M^'^K One may 
hope to infer information on the other entries from the other ^ processes. Many of 
them have been studied in the literature in this respect from both theoretical and 
experimental sides. Among them there are the decay — > /i+/i+7r^ |^, H, |], [1^, 
p!T| , the nuclear muon to positron |12| or to antimuon conversion, tri-muonium 
production in neutrino muon scattering and the process e'^p ul^l2X 

relevant for HERA [jl5[, as well as direct production of heavy Majorana neutrinos 



at various colliders fM. Unfortunately sensitivities of the current experiments 



searching for these processes are much less than in case of Oz//5/5-decay. The 



analysis made in the literature p, leads to the conclusion that if these processes 



are mediated by Majorana neutrino exchange then, except Oz//3/3-decay, they can 
hardly be observed experimentally. This analysis relies on the current neutrino 
oscillation data, and on certain assumptions related to the neutrino mass matrix. 
We will show that, despite the above conclusion being true for contributions of the 
neutrino states much lighter or much heavier than the typical energy of a certain 
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I/j process, there are still special windows in the neutrino sector which can be 
efficiently probed by searching for some of these processes. In the case of — >■ 
/i~^/x^7r~ decay this window lies in the neutrino mass range 245 MeV < niy. < 389 
MeV, where the s-channel neutrino contribution to the K"*" /i+/i+7r~ decay is 
resonantly enhanced, therefore making this decay very sensitive to neutrinos in 
this mass domain. If neutrinos with these masses exist, then from the present 
experimental data we can extract stringent limits on their mixing with z/^. 

Recently some phenomenological, cosmological and astrophysical issues of the 
intermediate mass neutrinos in the MeV mass region have been addressed [jl8| . 
This was stimulated by the attempts of explanation of the KARMEN anomaly 
in terms of these massive neutrino states |19|. Although recent data of the KAR- 
MEN collaboration have not con&med this anomaly, the possible existence of 
these massive neutrinos remains open, motivated by the idea of sterile neutrinos 
Ug required for the explanation of all the neutrino oscillation data including the 
LSND results. The sterile species may mix with the active ones i'e,^i,T to form 
massive states with a priori arbitrary masses. Their existence is the subject of 
experimental searches as well as cosmological and astrophysical constraints. 

The paper is organized as follows. In section 2 we discuss a model with 
sterile neutrinos and possible spectrum of massive neutrino states. Section 3 is 
devoted to some general features of constraints on the neutrino mass matrix, 
derivable from the ^ processes. In section 4 we give the theoretical framework 
for K"*" — > /i+zx+TT" decay, and then in section 5 discuss expected rates of this and 
other ^ processes in the light of the present neutrino observations. In section 6 
we study the possible contribution of massive neutrinos in the resonant domain 
of the yU+|U+7r^ decay and derive the constraints on the mixing of these 

neutrinos with u^^. Astrophysical and cosmological implications of these massive 
hundred MeV neutrinos are shortly addressed. 



2 Majorana neutrino mass matrix and neutrino 
counting experiments. 

Consider an extension of the SM with the three left-handed weak doublet neutri- 
nos u'l^i = (z/^g, z/^^, z/^^) and n species of the SM singlet right-handed neutrinos 
^Ri — i^Riy ■■■^Rn)- The general mass term for this set of fields can be written as 

- l^^Hz.- + H.c. = -l(z.' ( V + H.c. = 

2 2^ ^' \ Ml Mr J \ i^'^ J ^ 
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2 3+n 

~ -^^rn^ii^iUi + R.c. (1) 

i=l 

Here M.l,M.r are 3x3 and n x n symmetric Majorana mass matrices, Md 
is 3 X n Dirac type matrix. Rotating the neutrino mass matrix by the unitary 
transformation 

U^M^'^^U = Diag{m^i} (2) 

to the diagonal form we end up with n + 3 Majorana neutrinos Vi = Ul^v'^ with 
the masses mi,i. In special cases there may appear among them pairs with masses 
degenerate in absolute values. Each of these pairs can be collected into a Dirac 
neutrino field. This situation corresponds to conservation of certain lepton num- 
bers assigned to these Dirac fields. 

The considered generic model must contain at least three observable light 
neutrinos while the other states may be of arbitrary mass. In particular, they 
may include hundred McV neutrinos, which we will consider in section 6. Presence 
or absence of these neutrino states is a question for experimental searches. 

Let us point out that the presence of more than three light neutrinos are not 
excluded by the neutrino counting LEP experiments measuring the invisible Z- 
boson width Vinv Actually it counts not the number of light neutrinos but the 
number of active flavours. To see this let us write down the Zvu interaction term 
as 

n+3 n+3 

K^^a ^ UamU*^T7;;-f^iym = Y ^mn'^^T/.'^m, (3) 

a=e,n,T a=e,/i,Tm,n=l m,n=l 

where the last two expressions are written in the mass eigenstate basis. For the 
case of only three massive neutrinos one has Vmn — ^mn ^ ^ consequence of 
unitarity of In general Vmn is not a diagonal matrix and flavor changing 

neutral currents in the neutrino sector become possible at tree level. However if 
all the neutrinos are signiflcantly lighter than Z-boson with the mass Mz their 
contribution to the invisible the Z-boson width is 

n+3 

r,n.= Y l^mnl'rf^^rf^ Y Ma/? = 3r^^, (4) 

m,n=l a,P=e,ix,T 

where P^^ is the SM prediction for the partial Z-decay width to one pair of 
light neutrinos. This chain of equalities follows again from the unitarity of Uan- 
Thus, independently of the number of light neutrinos with masses m,^ << Mz/2 
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the factor 3 in the last step counts the number of weak doublet neutrinos. This 
conclusion is changed in the presence of heavy neutrinos N with masses Mn > 
Mz/2 which do not contribute to Fj^^,. In this case the unitarity condition is no 
longer valid and the factor 3 is changed to a smaller value. 

Having these arguments in mind we introduce in section 6 neutrino states 
with masses in the hundred MeV region. These states can be composed of sterile 
and active neutrino flavors as described in the present section. 

3 Constraints from fj -processes. General pat- 
tern. 

Let us examine some generic features of those constraints on the Majorana neu- 
trino mass matrix which can be derived from ^ -processes. 

Majorana neutrino mass terms in Eq. (||) violate lepton number conservation 
by two units AL = 2 and, thus, can induce -processes with AL = 2n. W-boson 
loops can be used to convert AL 7^ from the neutrino to the charged lepton 
sector. Therefore ^ -processes offer one of the most straightforward ways to test 
the Majorana nature of neutrinos and extract information on the neutrino mass 
matrix M^'^^. 

Note that at energies below the new physics thresholds only AL = 2n can 
be realized provided that baryon number is conserved AB = 0. This is a simple 
consequence of the Lorentz invariance and the fact that the spinor SM fields are 
represented only by leptons and quarks. To prevent AB ^ one has to contract 
the Lorentz indices of the external lepton fields only with each other without 
involving quark fields. Thus only ^ -processes with even number of external 
leptons, i.e. AL = 2n processes, can proceed at these energies. This means that 
any ^ -process in this energy domain is related to the Majorana neutrino mass 
receiving contributions from virtual Majorana neutrino exchange. Certainly only 
AL = 2 processes are of practical interest. The Majorana neutrino exchange 
contribution to the rate of a AL = 2 process with two external leptons lilj can 
be written schematically as 



The function in the absolute value brackets originates from the ^ Majorana neu- 
trino propagator (0|T(z/(a;)z/-^(?/))|0). Since the only source of ^ in the neutrino 
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sector is given by the Majorana neutrino masses rrij^k the decay rate (^) vanishes 
when m,yk = 0. The elhpsis in this equation denote terms whose exphcit form 
is irrelevant for the present general discussion. In Eq. G{z) is a smooth 
positively definite smearing function which depends on the particular process, sf 
are the limits of integration determined by the masses of the external particles 
involved in the process. The sign +(— ) in the denominator corresponds to the 
t-(s-) channel neutrino exchange. For the s-channel contribution the total neu- 
trino width r,^ has to be taken into account if masses in the resonant region of 
s-channel exchange Si < my < Si are considered. The non-zero neutrino decay 
widths Vyk can be introduced via the substitution m^fc rriyk — {i/'^)^uk- 

From the form of Fj^ one can infer that as a function on neutrino masses 
iriyk it has a maximal value V^"-^ for certain configuration of neutrino masses. 
This observation leads to the conclusion that the sensitivity F^^^ of a concrete 
experiment searching for ^ must satisfy the condition r^^^ < F™""^ otherwise 
no information on neutrino contribution is derivable. The experiment having 
passed this condition provides certain constraints. Assume that neutrino fields 
can be divided into light z/j and heavy Ni states with masses niyi « and 

\J~sl << Mtv respectively. Then in this "light-heavy" neutrino scenario the Eq. 
(PD can be approximately rewritten as 



m^AN ± Re 



moAyN (6) 



where the dimensionless coefficients Ai can be obtained for a concrete process 

from the formula such as Eq. (^. In the above equations niQ ~ -y/s^ is a typical 
scale of the ^ process under consideration. 

The average masses in Eq. (j^) are determined in the standard way 

1 > UikUjk 



k=light ^ k=heavy 



'^N k=heavy ^Nk 

Summation over light and heavy neutrinos implies masses rriyk << \[^ and 

A^Affc >> \[^ respectively. 

An experimental constraint on the rate of certain ^ process derived from its 
non-observation at the experimental sensitivity F^^^ can be translated with the 
aid of Eq. (|^) into the bounds: 
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However this is only possible if the experimental sensitivity satisfies the consis- 
tency conditions 



Exp{v) « y Si ^ rrio, Exp{N) ^ >> y sf ~ tjiq. (9) 
Otherwise experimental data can not be translated into the constraints (p|) as it 



was done, for instance, in Refs. [T^, |T5|. If the consistency conditions (y) are 



not satisfied one has to use the initial formula (|^). 

The following remark is in order. If all the neutrino states are light satisfying 
nT'uk « ^/sl then the following relation takes place 

(m,),, = M\f (10) 

This relation is not true if there are heavy neutrino states with masses not sat- 
isfying the condition m,y « \fs^. According to Eq. (|^) they do not contribute 
to {mu)ij measured in some ^ process. Therefore a concrete ^ process can give 
direct information on the entry A^lj^ of the Majorana neutrino mass matrix only 
under the assumption that all the neutrino masses are small compared to a typ- 
ical scale of this process m,yf^ « ttiq ~ or assuming that the heavy states 
are sterile. From this point of view ^ processes with larger typical scales mo are 
preferable. 

In the subsequent sections we will study concrete ^ processes from the view- 
point of their ability to probe neutrino properties. We will consider the conven- 
tional neutrino spectrum with three light neutrinos plus a number of kinematically 
unattainable heavy states as well as a model with additional intermediate mass 
neutrinos in the hundred MeV domain. 



4 K-meson neutrinoless double muon decay 

Here we shortly outline the theoretical framework for the — /i+/x+7r^ decay. 



Recently this ^ process attracted attention p], |Tl[| as a possible probe of the 
neutrino sector complimentary to other known processes. 

In the SM extension with Majorana neutrinos there are two lowest order 
diagrams, shown in Fig.l, which contribute to the K"*" — >■ fi'^fi~^7r~ decay. We 
concentrate on the s-channel neutrino exchange diagram in Fig. 1(a) which plays 
a central role in our analysis. The t-channel diagram in Fig. 1(b) requires in 
general a detailed hadronic structure calculation. In Ref. this diagram was 
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evaluated in the Bethe-Salpeter approach and shown to be an order of magnitude 
smaller than the diagram in Fig. 1(a), for light and intermediate mass neutrinos. 

The contribution from the factorizable s-channel diagram in Fig. 1(a) can be 
calculated in a straightforward way, without referring to any hadronic structure 
model. A final result for the K"*" fi^fi~^7c~' decay rate is given by |Tl[| 



"1 

r(K+ ^ /i+/i+7r-) = cj dsi 



E 



k Si - m^fc 



:iii 



K til Si-mlJ_ \ 



fm 



S2 - mf, 



The unitary mixing matrix Uij relates u'^ = UijUj weak u' and mass z/ neutrino 
eigenstates. The numerical constant in Eq. (|ll]) is 



c={G'^/32)in)-'fJlml\V^£\V. 



|2 

' US] ) 



(12) 



where = 1.28 = 0.668 and rriK = 494 MeV is the K-meson mass. 

The functions G{z) and H{zi,Z2) in Eq. (|ll]) after the phase space integration 
can be written in an explicit algebraic form 



h+-{z)h^^{z) -xlh^+{z) 



2 I / 2 \2 



(13) 



H{zi,Z2) = h^^{zi)h__{z2) + xl[r+{ziZ2) - xlt{zi,Z2, 1)] - r_{ziZ2)t{zi, Z2,X^). 



Here we defined xi = rrii/m^ and introduced the functions 

h±±iz) = z±xl±xl, r±{ziZ2) = Z1Z2 -xl± xj, 

t{zi, Z2, Z3) =zi + Z2- 2zl 0(z) = x% z)\^'\z, xj, xl) 

with A(x, y, z) = x"^ + y'^ + z^ — 2xy — 2yz — 2xz. 
The integration limits in Eq. ( pT]) are 



= m\{x^ + x^f 



m 



2y 



2y{l + xl)-{l + y- xl)h^+{y) ± 0(y) 



(14) 



(15) 



with y = Si/m?^. 
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Assuming that neutrinos can be separated into light Uk and heavy Nk states, 
with masses m,yi << \J~s^ and \J~^ « M^k we can rewrite Eq. ([TTD in the 
approximate form (|^) with the dimensionless coefficients 

A = 4.0 X 10"^\ An = 7.0 x 10"^^ A^n = 1-7 x IQ-^K (16) 

With these numbers we can estimate the current upper bound on the K"^ —>■ 
fi'^fi'^TT' decay rate from the experimental data on other processes. 



5 "Light-Heavy" neutrino scenario 

Here we assume that all the neutrino mass eigenstates can be divided into very 
light Ui and very heavy Ni states with masses, respectively, much smaller and 
much larger than the characteristic energy scale mo of the studied ^ process. 
Let us consider in this "Light-Heavy" neutrino scenario several typical exam- 
ples of ^ processes and estimate their ability to constraint the average masses 
{1/Mis[)ij as well as their possible rates. 
At present the highest experimental sensitivity to the Majorana neutrino con- 
tribution has been achieved in neutrinoless double beta decay (Oi//3/3) (A, Z) 
[A, Z + 2) + 2e~. A typical scale of this process is set by the nucleon Fermi 
momentum itlq ^ Pp ^ lOOMeV. The current constraints from OujSlS decay are 

i|] 

|(m,)ee| <0.2-0.6eV, | (-^)ee| < (9.0 • lO^GeV) "\ (17) 

The uncertainty relates to the uncertainties in the nuclear matrix elements and 
treatment of the background conditions. The ffist constraint in Eq. ([17D , as- 
suming that all the neutrinos are much lighter than 100 MeV, provides a direct 
constraint on the M.^^) entry of the Majorana neutrino mass matrix. Evidently 
these constraints satisfy the consistency conditions (^. 

Experiments searching for the other ^ processes have not yet reached enough 
sensitivity to establish meaningful constraints directly on the neutrino mass ma- 
trix elements. 

For instance, experiments on the muon to positron nuclear conversion /i~ + 
(A, Z) — > + (A, Z — 2) in ^'^Ti give at current sensitivity the following upper 



bound on the branching ratio |21 



R(/^- ^ e+) = ^^Jl^f ^ ^^\'^\ < 1.7 ■ 10-^^ (90%CL) (18) 
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Assuming that all the neutrinos are much lighter that the typical energy scale 
fno ~ "^fi = lOSMeV of this reaction one finds the bound 

IW^el < 17(80)MeV (19) 



for the proton pairs of the final nucleus in the singlet (triplet) state [|T2|. This 
constraints are marginal from the viewpoint of the consistency condition (^. 

Direct searches for K"*" — /i^/i+vr" decay by E865 experiment at BNL p3[ 
give 

T^.. = ^^^Y(^r^df^-^-^''^^~' (90^°CL). (20) 

Applying to this case the approximate formula @ one gets the limit 

Km,)^J < SOOGeV (21) 

which makes no sense because it does not satisfy the condition (H) with the typical 
energy scale of this process mg ~ = 494MeV. Thus the approximate formula 
(^ is not applicable to the present experimental situation for the K+ fi^fi^n^ 
searches. A similar picture holds for all known ^ processes searched for in various 
experiments (for other examples see Refs. ||17|| ). 



Viewing these processes from the stand point of neutrino observations one 
finds that except for the Oz//3/3-decay they have very small rates in the "Light- 
Heavy" neutrino scenario. 

Atmospheric and solar neutrino oscillation data demonstrate << (leV)"^ 
suggesting that all the neutrino mass eigenstates are approximately degenerate at 
the 1 eV scale . This observation in combination with the tritium beta decay 



endpoint allows one to set upper bounds on masses of all the three neutrinos ||2^ 



^e,^J.,T < 3eV. Thus in the three neutrino scenario one derives 

\{mu)ij\<9eV for i,j = e,fi,T. (22) 

This is much lower than the existing constraints on this quantity from ^ processes, 
except for Ou (3 jS-decaj which gives a significantly more stringent upper bound 
(plTl). With the constrain (|2^) we can predict the rate of various ]/j processes. 
Let us substitute the upper bound (p^ ) into the formula (|^) written for the 
—* fi'^fi^7T~ decay with coefficients given in Eq. (|T6|) . This gives rise to the 
following extremely small branching ratio 



7^^^ = ^ ^ ^ < 3.0 X 10"^° (3 light neutrino scenario). (23) 

r(K"'" — >■ all) 
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Assume there exist in addition heavy neutrinos N with the masses in the GeV 
region. Using the current LEP hmit on heavy unstable neutral leptons Mjv > 54.4 
GeV [0, we get 

\{M^')^,\<n{5AAGeVr\ (24) 

where n is the number of heavy neutrinos. 

This limit being substituted in Eq. together with the limit ( ]2^ ) results in 
the upper bound 

T^fifj. < 1-0 X 10^^^ (3 light + 1 heavy neutrino scenario). (25) 

Comparison of the theoretical predictions in Eqs. (|23|) , (^) with the experi- 
mental bound in Eq. (^) clearly shows that both cases are far from being ever 
detected. A similar conclusion is true for the other ^ processes except Ou^P 
decay. 

On the other hand experimental observation of these processes at larger rates 
would indicate some new physics beyond the SM, or, as we will see for the case 
of the K"*" — >■ /i+zi+TT" decay, the presence of an extra neutrino state uj with mass 
in the hundred MeV domain. As we discussed in section 2, the extra massive 
neutrino states Uj can appear as a result of mixing of the three active neutrinos 
with certain number of sterile neutrinos. These massive neutrinos are at present 
searched for in many experiments |^7|. The uj states would manifest themselves 
as peaks in differential rates of various processes, and can give rise to significant 
enhancement of the total rate if their masses lie in an appropriate region. 



6 Hundred MeV neutrinos in fi^fi^n -decay 

Assume there exists a massive Majorana neutrino Uj with the mass rrij 

^ 245 MeV < < ^ 389 MeV. (26) 

In this mass range the s-channel neutrino exchange diagram in Fig. 1(a) absolutely 
dominates over the t-channel diagram in Fig. 1(b), independently of hadronic 
structure. Here the diagram in Fig. 1(a) blows up because the integrand of the 
first term in Eq. (|TTp has a non-integrable singularity at s = m|. Therefore, in 
this resonant domain the total z/j-neutrino decay width Fj^j has to be taken into 
account. This can be done by the substitution rrij nij — {i/2)Vyj. 
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The total decay width T^j of the Majorana neutrino uj with mass in the 
resonant domain (^) receives contributions from the following decay modes: 

e~^7r~, e~7T~^, /^^vr", /U~7r+, 
— > <j e+e'u^, e+/i-i/;^, /i+e^z/^, fx^fi'u^^ , (27) 
e"e+z/e, e"/i+z/^, /i"e+z/e, fi'jj^^iyf,. 

Since i^,- = z/| it can decay in both z/j ^ l~X{AL = 0) and z/j ^ /"'"X^(AL = 2) 
channels. Calculating partial decay rates we obtain |ll|] 



r(u, In) = |f/,^f ^/>3F(y,,y.) ^ \Ui,M\ (28) 
1927r3' 



^ /i/2^) = |t/,^^f _^m5if(2/u, yi2) ^ \UiJ'Ti^'\ (29) 



where i/i = mi/rrij and 

i/(x, y) = 12 f-{z- y')il + x'- z)X'/\l, z, x2)A^/2(0, y^, ^). (31) 



2; 



21 



The integration limits are zi = yf^, 22 = (1 — y/J^ and F(0,0) = H{0,0) = 1. 
Summing up all the decay modes in (|^) one gets for the total uj width 

= 2|f/^,f (r^^) + r^'^') + r^'^^)) + 2|f/e,f (rf^ + r^^'^ + r^"^^) = 
= |f/^,fr('^) + |?7e,fr(^). (32) 



In the resonant domain (|26|) T^j reaches its maximum value at rrij = ysf. As- 



suming for the moment \Ufj_j\ = \Uej\ = 1, we estimate this maximum value to 
be Ti^j ^ 4.7 X 10~^° MeV. Since Tj^j is so small in the resonant domain (^61) the 
neutrino propagator in the first term of Eq. ([TT|) has a very sharp maximum at 
s = m|. The second term, being finite in the limit Ti^j = 0, can be neglected in 
the considered case. Thus, with a good precision we obtain from Eq. (|TT]) |Tl 



mj\Uf,j\ 



V 



r-(K+ /iV+TT-) ^ 071^(^0) rf^ TT (33) 

with zo = {'Trij/fT^KY ■ This equation allows one to derive, from the experimental 
bound of Eq. (PU[), constraints on the Vj neutrino mass mj and the mixing matrix 
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elements f/^j, Uej in a form of a 3-dimensional exclusion plot. However one may 
reasonably assume that \Uf^j\ ~ \Uej\- Then from the experimental bound (^) we 
derive a 2-dimensional rrij — If/^^jP exclusion plot given in Fig. 2. For comparison 
we also present in Fig. 2 the existing bounds taken from ||26|. As shown in 



the figure, the experimental data on the n^fi^n decay exclude a region 

unrestricted by the other experiments. The constraints can be summarized as 

< (5.6 ± 1) X 10"^ for 245 Mev < mj < 385 MeV, (34) 

The best limit < 4.6 x 10"^ is achieved at mj ^ 300 MeV. Note that these 

limits are compatible with our assumption that \Ufj_j\ ~ \Uej\ since in this mass 
domain, typically \Uej\'^ < 10"^ 



The constraints from //"''/i^vr" in Fig. 2 and Eq. can be signifi- 

cantly improved in the near future by the experiments E949 at BNL and E950 
at FNAL |2^. It is important to notice that in the resonant domain we have 



Yres^-^+ _^ ~ |f/|2, while outside r(K+ ^ fi+fi+n-) ~ \U\^. Thus in 

the resonant mass domain the K+ fi'^fi^n^ decay has a significantly better 
sensitivity to the neutrino mixing matrix element. In forthcoming experiments 
the upper bound on the ratio in Eq. (pOD can be improved by two orders of mag- 
nitude or even more. Then this experimental bound could be translated to the 
limit |t/u7p < 10^^^ and stronger. 



7 Hundred MeV neutrinos in astrophysics and 
cosmology 

It is well known that massive neutrinos may have important cosmological and 
astrophysical implications. They are expected to contribute to the mass density of 
the universe, participate in cosmic structure formation, big-bang nucleosynthesis, 
supernova explosions, imprint themselves in the cosmic microwave background 
etc. (for a review see, for instance, Ref. |29|). This implies certain constrains 



on the neutrino masses and mixings. Currently, for massive neutrinos in the 
mass region (|26|), the only available cosmological constraints arise from the mass 
density of the universe and cosmic structure formation. 

The contribution of stable massive neutrinos to the mass density of the uni- 
verse is described by the "Lee- Weinberg" VLyh'^ — irti, curve. From the requirement 
that the universe is not "overdosed" this leads to the two well know solutions 
iTLi, < 40eV and rrii, ^lOGeV which seem to exclude the domain Eq. (pBf). How- 
ever for the unstable neutrinos the situation is different. They may decay early to 
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light particles and, therefore, their total energy can be significantly "redshifted" 
down to the "overdosing" limit. Constraints on the neutrino life times t^^ and 
masses nij in this scenario were found in Ref. |^D[. In the mass region we 
have an order of magnitude estimate 

< (~ 10^'^)sec Mass Density limit (35) 

Decaying massive neutrinos may also have specific impact on the cosmic struc- 
ture formation introducing new stages in the evolution of the universe. After they 
decay into light relativistic particles the universe returns for a while from the 
matter to the radiation domination phase. This may change the resulting density 
fiuctuation spectrum since the primordial fiuctuations grow due to gravitation 
instability during the matter dominated stages. Comparison with observations 
leads to an upper bound on the neutrino life time [^|. In the mass region (|26|) 
one finds roughly 

T^. < (~ 10^) sec Structure Formation limit (36) 

On the other hand, on the basis of formula (p^), assuming ~ |f^ejP < 

4.6 X 10~^ as in Eq. (^4]), we find conservatively 

10~^sec < Ty., Theoretical limit. (37) 

Thus massive neutrinos with masses in the interval (|2^) are not yet excluded 
by the known cosmological constraints (^), (^) and there remains a wide open 
interval of allowed mixing matrix elements: 

(~ 10-i«) < |f/^,f , |f/e,f < (~ 10-^). (38) 

Big-bang nucleosynthesis and the SN 1987A neutrino signal may presumably 



lead to much more restrictive constraints |T^. Unfortunately, as yet the analysis 
[0 of these constraints does not involve the mass region (^). It may happen 
that these constraints, in combination with our constrains in Eq. (0), close the 
window for neutrinos with masses in the interval (PB]). Then the only physics left 
to be studied using the K+ — > /i+z^+yr" searches would be physics beyond the SM 
other than neutrino issues. Nevertheless, significant model dependence of all the 
cosmological constraints should be carefully considered before such a determining 
conclusion is finally drawn. 
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8 Conclusion 



We analyzed some generic properties of AL = 2 lepton-number violating pro- 
cesses and the constraints derivable from them on neutrino masses and mixing 
matrix elements. We discussed consistency conditions for experimental bounds 
when these bounds can be translated into the upper limits on the average neu- 
trino mass {mj) or the inverse average mass (I/Mtv). We found that excepting 
the neutrinoless double decay other AL = 2 processes are unable to provide us 
with sensible constraints on these quantities. Using the neutrino oscillation data, 
the tritium beta decay and the LEP searches for the heavy neutral lepton we esti- 
mated constraints on their rates in the scenario with three light and several heavy 
neutrinos. Typical values of these rates are far from being reached experimentally 
in the near future. 

We studied the potential of the K-meson neutrinoless double muon decay as 
a probe of Majorana neutrino masses and mixings. We found that this process is 
very sensitive to the hundred MeV neutrinos Vj in the resonant mass range (pGl). 
We analyzed the contribution of these neutrinos to the K"*" /i^/i+vr" decay rate 
and derived stringent upper limits on Majorana neutrino mixing matrix element 
|f/^jP from current experimental data. In Fig. 2 we presented these limits in the 
form of a 2-dimensional exclusion plot, and compared them with existing limits. 
The — > fi'^fi~^7i~ decay excludes a domain previously unrestricted experimen- 
tally. We stressed that the known astrophysical and cosmological constraints do 
not yet exclude hundred MeV neutrinos satisfying these K"*" — > /i+/i+7r~ con- 
straints. 

Finally, we notice that the decay K'^ fi^fi^n^ can in principle probe lepton 
number violating interactions beyond the standard model. However, according 
to recent studies supersymmetric interactions both with and without R- 

parity conservation seem to be beyond the reach of the K-decay experiments. 
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(a) (b) 

Figure 1: The lowest order diagrams contributing to K"*" — > ^'^fi~^TT~ decay. 
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Figure 2: Exclusion plots in the plane |f/^jp— mj. Here Ufj^j and rrij are the heavy 
neutrino i^j mixing matrix element to z/^ and its mass respectively. Domains above 
the curves are excluded by various experiments according to the recent update 
in Ref. [^. Region excluded by K"*" — > fj,'^fi~^7i~ decay ||Tl| covers the interval 



249MeV< mj < 385MeV and extends down to \U^j\'^ < 4.6 x 10"^ 
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